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We recently described a general solution to the phase matching problem that arises when one wishes to 
perform an arbitrary number of nonlinear optical processes in a single medium [Phys, Rev. Lett. 95, 133901 
(2005)]. Here we outline in detail the implementation of the solution for a one dimensional photonic quasicrystal 
which acts as a simultaneous frequency doubler for three independent optical beams. We confirm this solution 
experimentally using an electric field poled KTiOP04 crystal. In optimizing the device, we find — contrary 
to common practice — that simple duty cycles of 100% and 0% may yield the highest efficiencies, and show 
that our device is more efficient than a comparable device based on periodic quasi-phase-matching. © 2008 
Optical Society of America 
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1. Introduction 

Three-wave mixing is a nonUnear optical process that 
can take place within dielectric materials having a non- 
linear x^^^ couphng coefficient. Such processes are used 
for a variety of optical frequency conversion applications. 
Usually, due to dispersion, the three interacting beams 
do not propagate in phase and so efficient energy transfer 
between them is prevented.^ One of the common meth- 
ods to solve this problem, called quasi-phase-matching 
(QPM), is to periodically modulate the sign of the rel- 
evant component of the nonlinear dielectric tensor at a 
period corresponding to the phase mismatch.^' ^ This ap- 
proach is very successful, but unless one is extremely 
lucky it is limited to the phase matching of a single 
optical process. In recent years, the need to simultane- 
ously phase match several different processes arose in 
numerous applications such as creation of multiple radia- 
tion sources,^ of multi-colored solitons,^ of multi-partite 
entanglement sources,^ and for all-optical processing.^ 
This need was addressed by developing ad hoc general- 
izations for the quasi-phase-matching procedure, based 
either on periodic structures in one dimension^ (for non- 
collinear processes) and two dimensions®"^° or specific 
quasiperiodic structures in one^'^^^^^ and two^^ dimen- 
sions. In a recent publication [16, henceforth LAB] we 
explained how to solve the most general problem of mul- 
tiple phase matching — designing a device to phase match 
an arbitrary set of processes, defined by any given set of 
phase-mismatch values. The LAB solution is based on 
the general observation that the phase matching prob- 
lem is a consequence of momentum conservation, and 
that in crystalline matter, i.e. matter with long-range or- 
der, momentum conservation is replaced with crystal- 
momentum conservation. Thus, all that one needs to 
do is to design a nonlinear photonic crystal (NPC) — 
whether periodic or quasiperiodic — whose Fourier trans- 
form is peaked at all the required mismatch wave vec- 



tors. Here we present the first experimental realization 
of a device using this general solution: a onc-dimensional 
three- wave doubler. 

Note that other schemes for multiple harmonic gener- 
ation have been demonstrated before. ^^"^^ Nevertheless, 
we choose this relatively simple application of the LAB 
solution, as it allows us to provide a detailed pedagogical 
description of the approach. Other than demonstrating 
that the LAB solution indeed works, we wish to clarify 
all the steps in the design process, so that others will be 
able to implement it as well. We stress that the solution 
is general and is not limited to such simple applications. 

2. Simultaneous Phase Matching of Several In- 
teractions 

We consider second order nonlinear optical interactions 
in which three beams couple through the nonlinear sus- 
ceptibility x^^''. For a planar process in which two con- 
stant undepleted beams, Ei and E2, give rise to a third 
output beam E3, its integrated field amplitude is given 

by 

EsiAk) = r / g{r) exp(iAk • r)cfr, (1) 
J A 

where A is the interaction area, and F is a parameter 
depending on the amplitudes of the incoming waves, 
on the indices of refraction of all three waves, on the 
strength of the relevant component of the nonlinear sus- 
ceptibility tensor dij, and on the interaction width W. 
For example, for sum frequency generation in MKS, 
F ~ uj1dijEiE2/i(?k-iW ^ where c is the speed of light in 
vacuum. The function g{r) gives the spatial dependence 
of the relevant nonlinear coupling coefficient, and Ak is 
the phase mismatch vector of the interacting waves. For 
sum- frequency generation Ak would be ki -f k2 — ks. 

It is clear from Eq. (1) that the intensity of the out- 
put beam is proportional to the Fourier spectrum of the 
function g{v), evaluated at the mismatch vector Ak. 
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an NPC by fabricating whole strips normal to the tiles 
along the line. 

Before starting we wish to point out that in special 
cases, the D phase mismatch vectors Ak^-'-' may be in- 
tegrally dependent. This means that one can use fewer 
than D wave vectors to generate the NPC and still have 
Bragg peaks at all D points. It is then a matter of 
choice whether to use the full set of dependent vectors — 
although, as pointed out by LAB, it may be difficult 
in this case to control the intensities of the peaks — or 
to prefer a smaller set of independent vectors. Here we 
keep all mismatch vectors and treat them as if they were 
integrally independent. 



Fig. 1. (Color online) Illustration of the LAB solution for 
designing a one-dimensional NPC for multiple collincar 
optical processes, using the dual grid method, (a) The 
required mismatch vectors, (b) The dual grid, in which 
each family of lines is shown with a different color, (c) 
Tiling of the real-space line according to the order in 
which lines of different families appear in the dual grid, 
(d) Associating a given duty cycle with each tiling vec- 
tor. Positively-poled segment are shown in blue, and 
negatively-poled segments are shown in white. 



Thus, if we wish to simultaneously phase-match a set of 
D three-wave optical precesscs, characterized by phase 
mismatch vectors Ak^-'-', j = 1, . . . D, we should design 
the spatial structure of g{r) so that its Fourier spec- 
trum is peaked at all the D mismatch vectors. For a 
single process, the standard QPM solution^' ^ is to de- 
sign a one-dimensional NPC with a period of 27r/|AA;|, 
for which there is a first order Bragg peak in the spec- 
trum at Ak. The LAB solution shows how to design 
an appropriate nonlinear photonic crystal — whether pe- 
riodic or quasiperiodic — such that its spectrum contains 
Bragg peaks at any given set of D mismatch vectors. 
The approach that LAB adopt for this purpose is based 
on the so-called dual-grid method, originally developed 
by de Bruijn^* and later generalized^^"^^ to become one 
of the standard methods for creating tiling models of 
quasicrystals. 

The reader in encouraged to consult LAB^^ for a com- 
plete and rigorous treatment of the most general two- 
dimensional multiple phase-matching problem, which we 
do not repeat here. Instead, we give a detailed demon- 
stration of the LAB solution in one dimension, where all 
the optical processes arc chosen to be collincar. In this 
case the implementation of the dual-grid method for the 
design of the NPC, as well as the experimental fabrica- 
tion of the device, are relatively simple. Nevertheless, we 
can still design nontrivial and interesting devices, such 
as the three-wave doubler, implemented here. The ba- 
sic idea is to find a set of one-dimensional tiling vectors 
a^'^\ i = 1,...D, with which we can generate a one- 
dimensional tiling of the line, whereby a tile is simply 
an interval on the line. We then convert the tiling into 



3. Designing a One-Dimensional Three-Wave 
Doubler 

We wish to design a device that will simultaneously 
phase-match three collincar second-harmonic-gcncration 
processes, for three different wavelengths in the fiber 
telecom C-band: 1530nm, 1550nm, and 1570nm. We 
choose to use the nonlinear crystal KTiOP4 and oper- 
ate at a temperature of 100°C. At these conditions, the 
phase-mismatch values for the three processes are:^^'^'^ 
A/c(^) = .263//TO-1, Afc(2' = .256AiTO-\ and Afc'^' = 
.249/Ltm~^ respectively. Thus, we need to design an NPC 
whose Fourier spectrum contains three coUinear wave 
vectors with these dimensions, as shown schematically 
in Fig. la. 

In what follows we describe the design of such a struc- 
ture, as a particular example of the LAB solution for D 
collincar processes. Generalizing from D — 3 processes to 
an arbitrary number D of processes, follows directly by 
replacing all 3-component and 2-component vectors be- 
low by D-component and (-D — l)-component vectors, re- 
spectively. Generalizing to higher dimensional processes 
requires use of the full solution, as described by LAB. 

A. Finding the tiling vectors 

To calculate the corresponding three collincar tiling vec- 
tors, a^'^\ i = 1, . . .3, we first construct a single three- 
component vector ki = {Ak'-^\ Ak^^\ Ak'-'^^) from the 
three given mismatch values. This vector spans a one- 
dimensional subspace of an abstract three-dimensional 
vector space. Wc complete it to an orthogonal basis of 
the three-dimensional space by adding two orthogonal 
vectors q2 and qa. These are, of course, not unique, and 
we choose them to be q2 = (0.6483,-0.3421,-0.3331) 
and qa = (-0.3421, 0.6672, -0.3240). Wc use these three 
vectors as the columns of a 3 x 3 non-singular matrix, 



(2) 



whose rows K*^^^ , j = 1, . . . 3, span the three-dimensional 
vector space as well. We then find the three dual basis 
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Fig. 2. An optical microscope image of the demonstrated 
NPC. The prominent elements correspond to the a*^^^ 
tiling vector, poled with a 100% duty cycle. Their width 
is 8.5/iTO. The distances between these elements are 
quasiperiodically ordered along with the a^^^ and a^'^^ 
tiling vectors, whose widths are S.lfim and 7.9/iTO, re- 
spectively, and which are poled with a 0% duty cycle. 



vectors, denoted 




b'-P bf 

,(3) ,(3) 
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(3) 



by solving the three-dimensional orthogonality relations 



(4) 



Each row of the matrix (3) is a dual-basis vector of the 
form A'^' = (a*^^' , b'-')). The a'-'-' are the three required 
coUinear tiling vectors, whose values arc calculated to be 
a(i) = 8.3946/im, a^^) = 8.1652/im, and a^^) = 7.95/im. 
The 2-dimensional vectors b'^-'^ can be used, as explained 
by LAB, to analytically calculate the Fourier transform 
of the NPC in order to determine the expected efficien- 
cies for the different nonlinear processes. 

B. Constructing the tiling 

If we were now asked to generate all points at integral 
linear combinations of the three tiling vectors we would 
get the unwanted outcome of a dense filling of the real 
line. To avoid this situation we construct the dual grid, 
whose topology determines which of the integral linear 
combination of the tiling vectors are to be included in 
the one-dimensional tiling. 

The dual grid is constructed by associating with each 
mismatch vector Afc'-'^ j — 1,...3, an infinite family 
of parallel lines separated by a distance Lj ~ 2n/Ak'^^\ 
as illustrated in Fig. lb. The set of all families together 
constitutes the dual grid. We use the freedom to shift 
each family from the origin by an arbitrary value of fjLj, 
where < < 1, so that lines from different families 
never exactly coincide. Because the mismatch vectors are 
independent over the integers, such a shift produces a so- 
called gauge-transformation^"'^! which has no effect on 
the resulting NPC. 



The rest follows immediately, as the required order of 
the tiles in the real-space structure is given by the order 
in which lines of different families appear in the dual grid. 
This is illustrated in Fig. Ic. This is the sense in which 
the topology of the dual grid determines the real-space 
tiling in this trivial one-dimensional setting. The duality 
is a statement that each line in the grid is associated 
with a tile, or interval, in the tiling; and each interval 
in the grid with a vertex of the tiling. In our example, 
approximately 800 lines are required in each family to 
generate a 1cm long one-dimensional nonlinear photonic 
quasicrystal. 

C. Building a nonlinear photonic crystal from the tiling 

To create an actual nonlinear photonic quasicrystal we 
modulate the relevant component of the nonlinear sus- 
ceptibility tensor x^^-* according to the constructed tiling. 
Technology usually permits us to use a binary modula- 
tion of x^^-* so that the actual crystal can be represented 
by a normalized function g(r) = ±1. The simplest rep- 
resentation from a theoretical standpoint would be to 
attach a thin strip of value g(r) = 1 to every vertex 
of the tiling, while assigning the background a value of 
g{r) ~ —1. This is equivalent to a simple convolution of 
the strip with delta functions at the vertices of the tiling, 
and therefore gives the simplest analytical expression for 
the Fourier transform of the function g{r)}^ Neverthe- 
less, it does not necessarily produce an optimal NPC — 
one in which the strongest Bragg peaks are associated 
with the three mismatch vectors. In general, one can use 
a numerical procedure in order to optimize the required 
Bragg peaks. See for example the treatment by Norton 
and de Sterke.^^ Here we want to give a few quick-and- 
simple solutions — in addition to the thin strips — that are 
worth trying if one does not wish to deal with numerical 
optimization. 

A second option would be to use strips whose widths 
are equal to the tile vectors and simply to change the 
sign of g{r) from one strip to the next. In this way, ex- 
actly half the tiles will give strips of value g{r) = 1, and 
the other half strips of value g{r) = —1. The generated 
set of strips would be analogous to an antifcrromagnctic 
quasicrystal,^^' whose Fourier transform could also be 
calculated analytically. We have found that this option 
does not yield an optimal NPC for this application. 

A third option, and the one which we actually imple- 
mented, is again to use strips whose widths are equal 
to the tile vectors, but this time associate a so-called 
duty cycle with each tiling vector. This is done by divid- 
ing each tiling vector into two segments, and assigning 
a value of g{r) = 1 to one segment, and g{r) = —1 to 
the other. The duty cycle is the fraction of each strip 
with g{r) = 1. This general procedure is shown schemat- 
ically in Fig. Id. It gives us the ability to perform simple 
spectral shaping. By varying the three duty cycles, as- 
sociated with the three tiling vectors, we can engineer 
the magnitudes of the Fourier coefficients of the three 
required Bragg peaks. What we actually find — contrary 
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Fig. 3. (Color online) Spectral shaping. Each panel shows 
the magnitude of the Fourier transform for a 1cm long 
NPC made to phase match the three collinear processes, 
described in the text. In each panel one of the tiling vec- 
tors is given a duty cycle of 100%, denoted as 1, and 
the remaining two a duty cycle of 0%, or 0. Each panel 
also shows a piece of the corresponding real-space repre- 
sentation of the NPC, where the smallest element size is 
8^m. 
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Fig. 4. (Color online) Normalized pump wavelength re- 
sponse. Experiment and simulation results of second har- 
monic generation as a function of pump wavelength. This 
figure corresponds to panel (a) of Fig. 3. 

to common practice'^' ^^'^^ — is that the optimal NPC's 
are obtained when wc use duty cycles of cither 100% or 
0%. These are also easiest to fabricate in terms of the 
required resolution as we associate a value of g(r) = 1 
or g{r) = — 1 to tiles as a whole. 

In the experimental image of the NPC, shown in Fig. 2, 
the tiling vector a(i) is associated with strips of value 
g{r) = 1, or a duty cycle of 100%, while the other two 
have a value of g{r) ~ —1, or a duty cycle of 0%. In 
Fig. 3 we show numerical calculations of the magnitude 
of the Fourier transform of g(r), for the three possible 
assignments of a value of g{r) = 1 to one tiling vector, 
and a value of g{r) = —1 to the other two. One clearly 
sees that in all cases there are pronounced Bragg peaks 
exactly where we want them to be, but the distribution 
of intensities changes as we vary the tiling vector that 
is assigned a 100% duty cycle. The magnitudes of the 
Fourier coefficients are comparable to the 2/7r ~ 0.6366 
figure of merit, which is the magnitude of the first order 
Fourier component for a one-dimensional periodic NPC. 
In fact, because the efficiency is measured in terms of 
energy transfer it depends on the Fourier intensity and 
on the square of the interaction length. Thus, the real 
comparison should be with the efficiency per process of 
a sequence of three periodic NPC's, each of length L/3, 
which is proportional to (i • -l)^ = 0.045. For the homo- 
geneous nonlinear photonic quasicrystal that we imple- 
ment here the lowest process efficiency is proportional to 
(1-0.23)2 = 0.053 and the highest to (1-0.365)2 = 0.133, 
both better than the composite periodic structure. 

The ability to shape the spectrum using different rep- 
resentation of the function g{r) can be very useful for 
cascaded processes that to some extent are better per- 
formed in sequence. For example, in multi-harmonic gen- 
eration, where the outputs of some of the processes serve 
as inputs for others, it might be beneficial to give prece- 
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Table 1. Conversion efficiencies. Maximum conversion efficiencies are obtained for the indicated wavelengths. In the 
simulation, the maximum efficiency of the third process was shifted by +0.4nm. 



Wavelength 


Experimental efficiency 
[1/Watt] 


Theoretical efficiency 
[1/Watt] 


Simulated efficiency 
[1/Watt] 


1531nTO 


1.19e-3 


1.5e-3 


1.46e-3 


1551nTO 


3.37e-4 


6.13e-4 


6.83e-4 


1571nTO 


3.18C-4 


5.49e-4 


6.43e-4 



dence to the efficiencies of the initial processes at the 
beginning of the NPC and gradually, along the inter- 
action length, change the balance in favor of the latter 
processes. 

4. Experimental Results 

The NPC was fabricated by electric field poling^^ of 
KTiOP04, and is shown in Fig. 2. The spatial modula- 
tion was performed along 1cm with a 2 x 1mm? cross sec- 
tion. To test the NPC we used a pump beam from a tun- 
able external cavity diode laser, followed by an Erbium- 
doped fiber amplifier, and a fiber polarization controller 
whose purpose is to polarize the beam to be perpendic- 
ular to the plane of modulation of the NPC. The pump 
beam was then chopped at a frequency of IkHz and fo- 
cused to a waist of 20/im in the middle of the crystal. Its 
power was varied during the experiment in the range of 
4 — AOmW. The NPC was kept at a constant tempera- 
ture of 100°C. The input pump wavelength was varied 
in the range 1528nm — 1577nm and the resulting second 
harmonic power was measured using a calibrated Sili- 
con photo-diode, followed by a lock-in amplifier. The re- 
sults were compared with a simulation employing a split 
step Fourier method, in which a non-depleted Gaussian 
beam was used as a pump. Both results exhibit excellent 
agreement, as can be seen in Fig. 4, and show that this 
device is indeed most efficient for second harmonic gen- 
eration of the three designated pump wavelengths. 

Note that for focused Gaussian beams^° the actual 
peak efficiencies are shifted by Inm towards longer wave- 
lengths, as compared to the idealized design using plane 
wave interaction. However, with temperature tuning the 
results can be shifted in both directions. For the wave- 
length range used here, the shift rate is O.lnm per degree 
centigrade. The conversion efficiencies were also calcu- 
lated theoretically using the Boyd-Kleinman formalism 
for focused Gaussian beams, propagating within a ho- 
mogeneous nonlinear crystal, modified to include the 
relevant Fourier coefficient. The nonlinear conversion ef- 
ficiencies for the wavelengths for which conversion was 
maximal arc shown in Table 1. The small discrepancies 
between theory and simulation are attributed to the as- 
sumption, made in the theoretical calculation, that each 
processes is affected only by its relevant Fourier compo- 
nent. In the simulation the propagating waves experience 
small but still non-negligible contributions from close-by 
spectral components. 



5. Summary 

We present the first experimental realization of a nonlin- 
ear optical device employing the general solution of the 
phase matching problem, given by LAB. The demon- 
strated device is a one-dimensional three-wave doubler. 
We show that by simply using 100% or 0% duty cy- 
cles, each of the three wave doubling processes exhibits 
high efficiency. Moreover, the efficiencies arc all higher 
than for a device employing a periodic modulation with 
the same overall interaction length. We demonstrate the 
ability to perform spectral shaping of the response by 
changing the duty cycles, associated with the different 
tiles. This allows us to strengthen certain processes at the 
expense of others. It should be stressed once more that 
the LAB solution is general and can support any set of 
nonlinear processes, not only for one-dimensional prob- 
lems but also in two or three dimensions, without any 
symmetry restrictions. In addition, this general method 
can be applied to a broad range of problems from gener- 
ation of radiation sources through all-optical processing 
to generation of entangled photons in quantum optics 
applications. 
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